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Abstract. We propose a scheme for quantum lògic with neutral àtoms stored 
in an array of holographic dipole traps where the positions of the àtoms can be 
rearranged by using holographic optical tweezers. In particular, this allows for 
the transport of two àtoms to the same well where an external control ficld is 
used to perform gate operations via the molecular interaction between the àtoms. 
We show that optimal control techniques allow for the fast implementation of the 
gates with high fidelity. 
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1. Introduction 

In the search for a suitable system for quantum information processing, certain 
requirements have to be met [II . such as scalability of the physical system, the 
capability of initializing and reading out the qubits, and the possibility of having a set 
of universal lògic gates. Neutral àtoms are one of the most promising candidates for 
storing and processing quantum information. A qubit can be encoded in the infernal 
or motional state of an atom, and several qubits can be entangled using atom-light 
interactions or atom-atom interactions. Schemes for quantum gates for neutral àtoms 
have been theoretically proposed, that rely on dipole-dipole interactions |3 El EI EJ or 
controlled collisions |üj [7\ |SJ |Ü] . Such schemes can be implemented in optical lattices 
with a controlled filling factor, as shown in ref. |10j where multi-particle entanglement 
via controlled collisions was demonstrated. 

Presently a major challenge is to combine controlled collisions with the loading 
and the addressing of individually trapped àtoms. Recently techniques to confino 
single àtoms in micron-sized EU or larger dipole traps have been 

experimcntally demonstrated. A set of qubits can be obtained by creating an array 
of such dipole traps, each one storing a single atom |15|. Gate operations require the 
addrcssability of individual trapping sites and reconfigurability of the array. Arrays 
of dipole traps, each containing many àtoms, were obtained using either arrays of 
micro-lenses EvS] or holograms |17| . 

Actually, holographic techniques allow one to realize arrays of very small dipole 
traps |18| . which can trap single àtoms. Holographic optical tweezers use a computer 
designed diffractive optical element to split a single collimated beam into several 
beams, which are then focused by a high numerical aperture Iens into an array of 
tweezers. Recently holographic optical tweezers for individual Rubidium àtoms have 
been implemented by using computer-driven liquid crystal Spatial Light Modulators 
(SLM) |19|. The advantage of these systems is that the holograms corresponding to 
various arrays of traps can be designed, calculated and optimized on a computer. As 
a consequence, the trap array can be (slowly) controlled and reconfigured by writing 
thesc holograms on the SLM in real-time. 

Here we want to combine such an holographic array with a fast moving tweezer, 
in order to implement quantum gates based on a state selective collision between two 
àtoms, by using a Feschbach resonance. Optimizing the control of the àtoms motion 
is then of crucial importance, and is the subject of the present paper. 

2. Quantum register with holographic dipole traps 

The present approach for neutral àtoms quantum gates is related to several schemes 
which have been proposed for trapped ions and it uses a quantum register 

made of individual àtoms stored in an array of holographic dipole traps. The àtoms 
encoding the qubit will be stored in this register, which can be slowly reconfigured to 
move the àtoms around, but does not allow fast precise motion, which is required to 
implement a controlled collision between two àtoms. As a consequence, the register 
has to be combined with one (or several) fast tweezers, which can rapidly move an 
atom from one place to an other. There are then several options : either there is an 
atom in the moving tweezer, which can be entangled and disentangled with the àtoms 
in the register ("moving head" scheme, similar to the one proposed in ref. EOI)- O ne 
can also consider a configuration with two tweezers, which catch two àtoms in the 
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register and bring them to interaction. 

The fast tweezer (or tweezers) consist of a làser beam passing through an acousto- 
optical modulator (AOM), which allow to control simultaneously the deflection and 
the intensity of the beam with high accuracy. In the present paper, we will consider 
only two such tweezers, each containing one atom, and we will show that a quantum 
gate can be implemented with high fidelity by using optimal control techniques. The 
parameters of the calculations will be inspired by the experiment described in ref. 
jl·ll 1121 1131 I19j . but the scheme may work as well in a large range of parameter 
vàlues. Typically, the size of the beam waist for the tweezer will be less than a 
micron, resulting in oscillation freqüències of 130 kHz in the radial directions, and 
about 30 kHz in the axial direction. In addition, we will assume that a standing wave 
is added along the propagation axis. This has two important consequences : first, 
the axial oscillation frcqucncy is increased up to a valuc which is typically close or 
above the radial oscillation frequency; second, it will confine the two àtoms within the 
same "pancake" , therefore maximizing the non-linear phase shift acquired during a 
controlled cold collision. In the following, we will also assume that the two àtoms have 
been prepared in the ground state of the tweezer. Though this was not implemented 
yet, it can in principle be done, by using either side band cooling, or evaporative 
cooling down to the single atom level |13| . 



3. Atom transport in a time-dependent double-well potential 

The transport mechanism is discussed in [22] for àtoms in a time dependent, optical 
súper lattice which has the form of a periòdic array of double well potentials. Here, 
however, we consider a one-dimensional system with a single double well potential of 
the form 

V(x, t) = -A{t) e-* 2 /2- 2 _ B{t ) e -(x+d(t)f/2 W 2 (1) 

In the geometry described above, it is sumcient to consider the one-dimensional case, 
where the position coordinate x corresponds to the distance between the two tweezers. 
We will thus assume that the motional state along the two other axis does not change 
during the transport process to be described (this point is further discussed later in 
this section). 

The location and the depth of the minima of the potential is dctcrmined by 
the time dependent control parameters A(t),B(t) and d(t). The time evolution of the 
motional degrees of freedom of a single particle in the trap is governed by the time 
dependent Schròdinger equation 

ihj t tP(x,t) = H(t)i>(x,t) (2) 

with 

H(t) = ~^^ + V(x lt ). (3) 
2m ax 

In the following, distances are measured in units of a harmònic oscillator length a = 
y/Ti/mui and energies in units of e = h 2 /2ma 2 . In case of 87 Rb and lo = 2ir x 100kHz 
this defines a length scale of a = 34 nm and an energy scale of e = 2-kTi x 50 kHz. 

As in E3, we assume that there is initially one atom in the ground state of 
each well and that the barrier is sufhciently high to prevent tunnelling between the 
wells. This allows to raise rapidly the left potential well, such that at time t = 
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Figure 1. The double well potential at (a) the initial time t = 0, (b) at an 
intermediate time í = T/2 and (c) at the end of the transport process í = T. The 

position is given in units of a = y/ h/mu> and the energy in units of e = fi 2 /2mn 2 
as described in the text. The horizontal, dashed lines indicate the eigenenergies 
of the system. The sòlid lines are the (real) eigenfunctions of H(t) corresponding 
to the moving and the register atom, i.e. 1/12(2:, t) and Vo(^,í), respectivcly. 
The shown potentials correspond to experimental parameters as described in 
(ÏT1IT21IÏ3) . 

the situation depicted in figure H^a) can be created: The lowest motional state of 
the left atom corresponds to the second excited state of the double well system while 
the right atom is in the ground state. The moving atom can then be adiabatically 
transported to the right well by lowering the left well and the barrier simultaneously 
(see figure Elb)), and raising the barrier again while the left well is further lowered 
leading to the final configuration at t = T shown in figure ^c) . During the wholc 
process the moving atom stays always in the second excited instantaneous eigenstate 
of the system which corresponds eventually to the first excited state of the right well 
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Figure 2. (a) The lowest instantaneous eigenenergies of the doublc-well system 
during the transport process (the ground state energy is set to zero). The 
transported atom corresponds to the sòlid line and the static atom to the zero 
line. (b)(c) The pulse functions A(t), B(t) and d(t) control the depth and the 
distance of the two wells. In (a)-(c), the parameters for í = 0, t = T/2 and t = T 
are the same as in figure HI 



while the register atom remains in the ground state. 

The adiabatic transport is possible since by choosing appropriate pulse functions 
A(t), B(t) and d(t), level crossings of the eigenenergies of the system during the process 
can be avoided. An example for such pulse functions, as well as the corresponding 
instantaneous eigenenergies of H(t), are shown in figure [21 In figure [2[a) the time 
dependent energy of the moving atom is given by the bold line (the ground state 
energy is set to zero). Back transport of the moving atom to its original position is 
obtained by time inversion of the pulses. 

In order to study the dynamics of the transport process we introduce the 
occupation probabilities 

2 



dxtp (x,t)ip n (x,t) 



(4) 



where ip n (x, t) with n = 0, 1, 2, ... is the nth instantaneous eigenfunction of the double 
well potential. The superscript A £ {M, R} indicates the wavefunction of the atom 
to be transported (moving atom) and the atom which is supposed to stay located at 
its well (register atom), i.e. ip M (x,t) [ip R (x,t)] is the solution of the time dependent 
single particle Schròdinger equation J5J) with initial condition ^jj ÀI (x,0) = -02(a;,O) 
[ip R (x,0) = ïpo(x,0)] as shown in figure^a). The fidelities of the processes are then 
given by F M = F 2 M (T) and F R = P H (T). 

For the example shown in Fig.|21we get F M = 99.7% for propagating the moving 
atom wavefunction and F R — 99.9% for propagating the register atom wavefunction 
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Figure 3. (a) Fidclity F M (t) corresponding to the moving atom during the 
transport for T = 250ft/e (i), T = 350h/e (ii), T = 400h/e (iii), T = 5007i/e 
(iv), T = lOOOfi/e (v), T = 5000fi/e (vi), (b) Fidelity F R (t) corresponding to the 
register atom during the transport for T = 500?ï/£. (c) Probabilities P^ ' 1 (t) (sòlid 
line), Pf(t) (dash-dotted line), P 3 M (í) (dashed line) and P 4 M (í) (dotted line) of 
finding the moving atom in the respective eigenstates during its transport. 



from í = 0to£ = T = 500?i/e. In the case of Rubidium this would correspond to a 
time T — 1.6 ms. 

Figure[3]shows the probability P 2 M if) that the moving atom remains in the second 
instantaneous eigenstate during the transport for various operation times T. For 
T = 500h/e the fidelity is always greater than 98.7% and the corresponding occupation 
probability P^(t) of the register atom (shown in figure E£b)) is always larger than 
99.9%. The corresponding probabilities of finding the moving atom in the ground, 
the first excited, the third excited and the fourth excited instantaneous eigenstates 
are displayed in figure E^c). The occupation probabilities of highcr excited states are 
smaller than 4 x 10 -5 and are not shown. An excitation energy of 100kHz, which 
would be required for the excitation of radial motional states, would correspond at 
least roughly to the eighth excited state along the axis of motion. The occupation 
probability of this (and higher) state(s) is found to be smaller than 5 x 10 -9 which 
justifies the one-dimensional model used in this paper. 

In the following we will discuss the influence of experimental imperfections, 
especially variations in the làser intensities, which are proportional to the pulse 
functions A(t) and B(t), and variations in the distance of the làsers creating the double 
well potential, which affect the pulse function d(t). Motivated by the experimental 
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conditions we assume variations of the pulse functions of the form 

d(t) = d(t) + Sdsm(Qt) (5) 

À(t) = A(t) +SAsm(flt) (6) 

B(t) = B(t) + SB sin(Oí) (7) 

with íl = 27r/lms. Assuming a drift Sd = lnm while SA = SB = leads for 



T = 500h/e to a slight reduction of the fidelity to F M = 99.4% and a higher fidelity 
F R . A variation of SA = O.le while SB = results in F M = 99.5% and F R = 99.9%. 

If A(t) and B(t) undergo the same perturbation, i.e. SA = SB, the shape of 
the potential Q does not change significantly if the variation is not too large (except 
for an approximately constant shift of the potential). Therefore, the level spacing as 
shown in figure^a) remains roughly the same and it is expected that the transport 
can be done as fast as without fluctuations. This behavior is confirmed by numcrical 
simulations: Assuming SA = SB = 10è, which corresponds to variations of the làser 
intcnsity of approximately 1%, we get the fidelities F M = 99.7% and F R = 99.9%. 

This analysis shows that the current transport scheme is relatively insensitive to 
noise which affects both parameters, A(t) and Bit), in the same way, while it is more 
sensitive to different perturbations in these parameters. In this case, level crossings 
in the energy diagram[3]can appear, leading to significant leakage into higher excited 
states, which would require a more sophisticated engineering to be controlled and will 
be a subject of future investigations. 

4. Quantum gates by optimal control of molecular interactions 

Performing gate operations requires a strong molecular interaction between àtoms. 
They can be coupled to molecular states either by means of Feshbach resonances [23! 
or through Raman photo-association làser pulses [21] • For the sake of concreteness, 
we focus here on Feshbach resonances - however, all of our arguments can be adapted, 
e.g., to Raman photo-association. We consider 87 Rb àtoms. 

Feshbach resonances occur when a bound molecular state |n) crosses the 
dissociation threshold for a state having the same quantum numbers \2'A\ while 
changing an external magnètic field B. Close to resonance, the scattering length 
varies as 

A{B) = A bg (l-^^j, (8) 

where A bg is a non resonant background scattering length, Bq is the resonant magnètic 
field, and A„ is the width of the resonance. The resonance energy varies almost lincarly 
with the field 

e n (B) = s n (B - B ), (9) 

with a slope s n . We are interested in the dynamics of such a system in a confmed 
geometry. Following [21], we shall model it by the effective Hamiltonian 

(B)\n)(n\ + ^2(vhiy\v)(v\ + V v \v)(n\ + h.c), (10) 

V 

where the |u)'s are the trapped relative-motion atòmic eigenstates of an isotropic 
harmònic oscillator trap having frequency v. The couplings to the resonance are 

V v = 2fivJ V4u +"3 a bg 5 n /ir (11) 
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Figure 4. Two-qubit gate operation via optimal magnètic field control: 
optimized field time dependence (top left); overlap between initial and evolved 
state (bottom left); accumulated two-particle phase tp (top right); decrease of the 
infidelity with increasing iterations (bottom right). 



with cibg = Ai,g S n = A n s n / (Tiu). In a different geometry, for instance in 
an elongated trap characterized by a ratio 7 between the ground level spacings in 
the transverse and in the longitudinal potential, the couplings can be calculated by 
projection on the corresponding eigenstates [221 • Accurate vàlues for the resonance 
parameters A„ and Bq, as well as for Af, g , are now available from both theoretical 
calculations and recent measurements 26 . 

The possibility of controlling the resonance energy via an external magnètic field, 
as described by Eq. JÜJ), provides a straightforward way to steer the interaction between 
the àtoms. Indeed, the coupling to a specific resonant state \n) is only effective for 
a particular entrance channel, i. e. a specific combinations of atòmic hyperfine states 
(that is, of logical qubit states in our case), while in general all other channels will 
be unaffected by the resonance. Thus the resonance-induced energy shift will cause 
a two-particle phase to appear only for that particular two-qubit computational basis 
state. 

We will identify our qubit logical states with the clock-transition states 

|0) = \F= l,m F = 0), 1 1) = \F = 2, mp = 0) . (12) 

The main advantage of this choice is that the qubit states are not sensitive to the 
magnètic field, and hence not subject to decoherence due to its fluctuations. We will 
use the resonance for the channel |00) oceurring around Bo = 685 G, having a width 
A„ = 16 mG. For obtaining a two-qubit gate, the magnètic field is ramped across Bq, 
and eventually tuned out of the Feshbach resonance again, getting the following truth 
table for the operation: 

|00) -> e iv |00), 
|01) - |01>, 
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|10)-|1Q), 

|11>-|11>> (13) 
where we included the phase <p accumulated by state |00) during the ramping process 
due to the interaction energy shift, whose value can be adjusted by controlling the 
magnètic field. If (p = ir, a C-phase gate between the two àtoms is obtained. Note 
that làser addressing of single qubits is never required throughout the procedure. 

The magnètic ramping process can be even performed non-adiabatically, provided 
that all population is finally returned to the trapped atòmic ground state. This can 
be accomplished via a quantum optimal control technique in analogy with the above 
discussion for the transport process. The control parameter in this case is the external 
magnètic field B. Care has to be taken in optimizing not only the absolute value 
of the overlap of the final state onto the goal state, but also its phase <p. Fig. 0] 
shows the optimization results for a trap with transverse freqüències of f 00 kHz and 
a longitudinal frequency of about 25 kHz, corresponding to the right well of Fig. ^ 
The final infidelity is about 2 x 10 -5 in this case. 

5. Outlook 

We have described a scheme using moving twcezcrs and statc-dcpcndcnt controllcd 
collisions, which is able to implement a quantum gate between two individual àtoms 
with a high fidelity. The sensitivity of the scheme to intensity or position fluctuations 
has been examined, and the controlled motion is found to be very tolerant to "common- 
mode" noise between the two tweezers. It is even relatively tolerant to differential 
noise, because the overall process is close to adiabatic, and designed in such a way as 
to avoid unwanted level crossings. 

The magnètic field which is used here to obtain the Feshbach resonance would 
be ultimately very advantageously replaced by an optical field [23 HE! > which can 
be switched on and off with high speed and precision, and which will not perturb 
the neighboring àtoms stored in the holographic array. Though the overall scheme is 
clearly not easy to implement, optimal control techniques as used here certainly help 
to make it closer to realistic. 
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